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ABSTRACT 

Rotating white dwarfs undergoing quasi-radial oscillations can emit gravita- 
tional radiation in a frequency range from 0.1 — 0.3 Hz. Assuming that the en- 
ergy source for the gravitational radiation comes from the oblateness of the white 
dwarf induced by the rotation, the strain amplitude is found to be ~ 10~^^ for a 
white dwarf at ~ 50 pc. The galactic population of these sources is estimated to 
be ~ 10'', and may produce a confusion limited foreground for proposed advanced 
detectors in the frequency band between space-based and ground-based interfer- 
ometers. Nearby oscillating white dwarfs may provide a clear enough signal to 
investigate white dwarf interiors through gravitational wave asteroseismology. 

Subject headings: gravitational waves — stars: oscillations — white dwarfs 



1. Introduction 



There are a number of gravitational radiation detectors, planned, under construction, 
and operational covering a wide frequency spectrum from ~ 10^^ Hz all the way up to 
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~ 10^ Hz. Expected sources of gravitational radiation include numerous astrophysical sources 
such as compact binaries, supermassivc black holes, and binary coalescence. In addition, 
there is an expected cosmological background of gravitational radiation arising from the 
very earliest times of the universe. The coverage of the spectrum is not complete and the 
gap between space-based interferometers (such as LISA) and ground-based interferometers 
(such as LIGO and VIRGO) has been proposed as a possible "clean window", devoid of 
continuous foreground sources, through which the cosmological background of gravitational 
radiation could be seen (Seto, Kawamura, & Nakamura 2001). Most binary white dwarf 
systems will coalesce before their gravitational wave frequency rises above 0.1 Hz, while more 
massive binaries such as double neutron star or black hole binaries will be sweeping through 
this band on the way to their eventual coalescence in the ground-based frequency band. 
However, rotating white dwarfs undergoing quasi-radial oscillations will emit gravitational 
radiation in this frequency band. These sources will be essentially monochromatic and long- 
lived. Given the number of white dwarfs in the galaxy, it is quite possible that this population 
will produce a confusion-limited foreground of sources in this frequency band that will mask 
the cosmological background. We propose a possible energy source for these oscillations and 
attempt to estimate the average signal strength. 



Quasi-radial oscillations of rotating white dwarfs were investigated in the early 1970's (Pa- 
poyan, Sedrakian, & Chubaryan 1972; Haroutyunyan, Sedrakian, & Chubaryan 1972) where 
the frequency spectrum of the fundamental oscillation mode for maximally rotating white 
dwarfs was determined. These stars are oblate due to their rotation and consequently 
they have a non-zero quadrupole moment. The oscillations add a time dependence to the 
quadrupole moment (Vartanian & Hajian 1977). The oscillation is described by assigning 
each mass element a time dependent coordinate given by Xa = a:°(l + r)smu}t) where rj 
and a constant. Thus, the reduced quadrupole moment is given by: 



where Q^p are the components of the quadrupole moment of the rotating oblate white dwarf 
in equilibrium and we have neglected terms of order 77^. Taking the axis of rotation to lie 
along the z-axis, the non-zero components of the quadrupole moment obey: 



2. Quasi-radial Oscillations 




g" (l + 277sina;i) 



(1) 



|0 



(2) 
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The power emitted in gravitational radiation is given by: 



5c5 



(3) 



and consequently one obtains: 



J = ^r^V |g°|'cos'a;i' = Jocos'u;t' (4) 

where the retarded time is = f — r/c for a source at distance r. 

To determine the wave form and the angular distribution of the radiation, we rotate to 
coordinates in which the wave vector lies along the 2;-axis and use the transverse-traceless 
gauge. Consequently, 

hx = Ky = 0, (6) 



where 6 is the angle between the wave vector and the axis of rotation of the white dwarf. 
We can express the strain amplitude in terms of the power by combining Eq. 4 with Eq. 5 
to obtain: 

/15GJ0 1 . 2/3 • frj. 

= \ sm t^smo;!. (7) 

V 2c-^ ruj 

If an energy source can be found to drive the pulsations, the rate at which power is put into 
the vibrations can be combined with the lifetime of the energy source to estimate the strain 
amplitude from an individual white dwarf and thus from the galactic population as a whole. 
We discuss a possible mechanism in the next section. 



3. Energy Source 

If there is no permanent source of energy to feed the gravitational radiation, the os- 
cillation energy will quickly radiate away in about 1000 years (Vartanian & Hajian 1977). 
Since the ultimate source of the gravitational radiation from the white dwarf is the oblate- 
ness arising from the rotation, we propose that the deformation energy of the white dwarf 
provides the energy to drive the oscillations. In this scenario, as the white dwarf spins down, 
it will transition from oblate to spherical. This transition will trigger starquakes which will 
feed the oscillations which drive the gravitational radiation. We assume that a part of the 
deformation energy will be converted into gravitational radiation while the remaining part 



-4- 



will dissipate through thermal, electromagnetic, and other channels. We use the technique 
of Sahakian, Sedrakian, & Chubarian (1972) to calculate the deformation energy. 

Prom numerical results, the dependence of the mass of the rotating and non-rotating 
configurations is found to be linear with the baryon number, so M = kN and Mq — koN, 
where the subscript "0" indicates the non-rotating configuration. The mass difference be- 
tween the rotating and non-rotating configurations with the same baryon number is a result 
of the additional energy of rotation [Wr{fl) — ^I^"^) as well as the potential energy due to 
crustal deformation {Wg{fl)), thus: 

Wg{n) = AMc^ - Wr(n), (8) 

where Q is the angular velocity of the white dwarf, / is its moment of inertia, and the mass 
difference in grams is given by 

AM ^{k- ko)N = 8.96 x IQ-^^N. (9) 

The appropriate parameters for maximally rotating white dwarfs were calculated in Haroutyu- 
nian, Sedrakian, & Chubaryan (1971) and Sahakian, Sedrakian, & Chubarian (1972). These 
parameters are presented in Table 1. To obtain the deformation energy for rotation rates less 
than fimax) "we use the fact that the results of Haroutyunian, Sedrakian, & Chubaryan (1971) 
were obtained using a linear expansion in the small dimensionless parameter fl"^ /SnGpc, 
where pc is the non-rotating central density. Consequently, we can write: 

Wg{n) = (-^) Wg(n^^). (10) 

\ ^ 'max / 

It now remains to determine the time scale, r, for a spin-down mechanism so that we can 
relate the power in gravitational radiation to the decrease in deformation energy by: 

Jo^P-^, (11) 
r 

where /3 ^ 1 is a "branching ratio" that quantifies the fraction of deformation energy which 
goes into gravitational radiation from the fundamental mode. 

One possible mechanism for the spin-down of a rotating white dwarf is magneto-dipole 
radiation torque, which occurs if the magnetic field is oblique (Schmidt, et al. 2001). Ob- 
servational data for 65 isolated white dwarfs indicates the magnetic field strength on the 
surface of these stars lies in the range ~ 3 x 10^ to ~ 10^ G (Wickramasinghe & Ferrario 
2000). If we define a to be the angle between the magnetic and rotation axes, the spin-down 
rate of the white dwarf is given by: 

n = -|yi-sin^a, (12) 
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where /i — BE? is the magnetic moment, B is the magnetic field strength, and R is the 
radius of the white dwarf. The characteristic time scale is then: 



T 



(13) 



We neglect the influence of gravitational radiation on the spin evolution because the ratio 
of the energy lost through gravitational radiation to that of the energy lost through the 
magnetic-dipole spin down mechanism is about 3.6/3. With this information, we can now 
estimate the gravitational radiation luminosity and strain amplitude for several white dwarfs. 



4. Strain Amplitudes and Luminosities 

The rotation rates for white dwarfs arc difficult to measure since these objects have little 
or no surface blemishes and gravitational broadening of their spectral lines overwhelms the 
expected rotational broadening (Kawaler 2003). Fortunately, some magnetic white dwarfs 
show time variability in their magnetic features which allows for their rotation rates to 
be inferred (Wickramasinghe & Ferrario 2000; Kawaler 2003). Magnetic white dwarfs are 
thought to make up roughly 1% of the white dwarf population and isolated magnetic white 
dwarfs tend to have an average mass (> O.95M0) that is well above that of nonmagnetic 
white dwarfs (~ O.STM©). There is evidence that the population of magnetic white dwarfs 
is bimodal with high mass stars (near the Chandrasekhar limit) having rotational periods 
between ~ 700 s and several hours (Wickramasinghe & Ferrario 2000). This may be evidence 
that such objects are the result of double degenerate mergers (Wickramasinghe & Ferrario 
2000; Dupuis, Vennes, & Chayer 2002), as has been proposed for EUVE J0317-855 (Ferrario, 
et al. 1997). 

The population of merged magnetic white dwarfs may be the most promising source 
of gravitational radiation from the mechanism of quasi-radial oscillations. To estimate the 
expected strain amplitude from both individual sources and the population as a whole, we 
identify eight candidates from the 65 known isolated magnetic white dwarfs from Wickramas- 
inghe & Ferrario (2000) as those with rotational periods less than one day. The properties of 
the candidates are listed in Table 2. We estimate the local space density of merged magnetic 
white dwarfs by taking the local space density of white dwarfs to be 0.003 pc~^ (Liebert, 
Dahn, & Monet 1988), the fraction of isolated magnetic white dwarfs to be 5.1% (Wick- 
ramasinghe & Ferrario 2000), and the fraction of isolated magnetic white dwarfs that are 
remnants of mergers to be 8/65. This results in a local space density of ps = 1-9 x 10~^ pc~^. 
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The expected distance to the nearest source is then found by: 

1/3 



47rp, 



~ 46 pc (14) 



Assuming a galactic distribution of white dwarfs to follow the disk population, we assign a 
density distribution of: 

p = (15) 

in galacto-centric cylindrical coordinates, with R = 2.5 kpc and h = 200 pc. Taking the 
solar location as = 8.5kpc and Zg = 0, we obtain pq = 5.5 x 10~^ pc~3 and a total number 
of = 8.6 X 10^ in the galaxy. 

Expected strain amplitudes, r, and r] for the eight isolated magnetic white dwarfs are 
given in Table 3. These are calculated by first determining r from Eqs. 12 and 13, where we 
have chosen sin^ a = 1/2. This is a reasonable assumption since a distribution in a that is 
cither uniform or spiked at a; = 7r/4 arc supported by observation (Schmidt & Norsworthy 
1991). Choosing a value for f3 with which to calculate the luminosity is somewhat problem- 
atic due to the dearth of observational data for pulsations in this frequency range. Most 
observations of white dwarf pulsations are in the mHz range, and consequently any pulsa- 
tions in the dHz range will be averaged out in an observation. We choose an upper bound to 
/3 by requiring that the largest Doppler broadening of spectral lines due to pulsations be less 
than thermal Doppler broadening. The most stringent constraint comes from PG 1031-1-234 
and yields f3 = 10^^. We note that the resulting pulsational Doppler broadening in the 
remaining seven white dwarfs is at least an order of magnitude below the thermal Doppler 
broadening. Thus, the luminosity is calculated from Eq. 11 and used with Eqs. 14 and 7 
to determine the strain amplitude where we have averaged over all orientations and taken 
P — 10""^. We have used the average mass of 0.95 Mq whenever the mass was undetermined 
and the average distance of 46 pc is the distance is unknown. The expected energy flux on 
earth (in erg s~^ cm~^) for a population made entirely of each type of white dwarf is also 
shown in Table 3. The flux, F, is calculated using 

F = ^^P^^f(^^\ (16) 
47r(3xl0i8)' ^ ^ 

where h = 200 pc and f{h) = 6.15 is calculated using Appendix A of Hils, Bender, & 
Webbink (1990). Finally, we note that a simple average of the strain amplitudes in Table 3 
gives h+ — 5.9 x 10~^^ and an average flux of F = 9.21 x 10"-^^ erg s"^ cm~^. The flux is 
spread out over a frequency band of ui — 0.12 to 1/2 — 0.32 Hz, and we can estimate an 
average strain amplitude for the galactic population of pulsating white dwarfs by using the 
angle and polarization averaged expression of Douglas & Braginsky (1979) and averaging 
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over the frequency range Ai/ — 1^2 — t^i to obtain: 



h. 



■+ave — 



Inz/a/z/i lAGF 

Av V TTC^ ' 



(17) 



which gives /i+ave 



= 8.35 X 10-2^ 



5. 



Conclusions 



We have shown that the galactic population of white dwarfs can produce a background 
of gravitational radiation in the frequency range of 0.12 — 0.32 Hz through quasi-radial pul- 
sations. The source of energy to drive these pulsations is found in the deformation energy of 
the white dwarf due to its rotation. This energy can be extracted from the white dwarf as 
it spins down. We have proposed that a population of isolated magnetic white dwarfs which 
are the remnants of merged double degenerate binaries can be sources of this gravitational 
radiation. Estimates of the signal strength over the frequency band of interest indicate that 
this population may be comparable in strength to the level of the stochastic cosmological 
background of gravitational radiation predicted by standard inflationary models. In addition, 
some white dwarfs may be near enough that their signals will stand above the background 
and their individual parameters can be measured. If the frequency and amplitude of the 
gravitational radiation from one of these nearby single white dwarfs can be measured, then 
it will provide an opportunity to test models of white dwarf interiors and give a unique op- 
portunity to check the main parameters for white dwarfs and to understand their dissipation 
mechanisms. 
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Table 1. Structural Parameters of Maximally Rotating White Dwarfs 



Pc(6) 


iV(57) 


M/Mq 


^max 


/(48) 


^(48) 


R{8) 




W^r(49) 


^9(49) 


u 


2.403 


0.4997 


0.5945 


0.196 


128.0 


20.48 


10.93 


4.799 


0.246 


4.55 


0.758 


19.38 


0.8398 


0.9993 


0.476 


88.6 


14.27 


7.342 


8.066 


1.00 


7.06 


0.794 


157.7 


1.0695 


1.2731 


1.063 


39.5 


4.766 


4.625 


10.27 


2.23 


8.05 


1.513 


866.1 


1.1340 


1.3502 


2.042 


15.9 


1.554 


3.044 


10.90 


3.32 


7.58 


1.985 


2586 


1.1261 


1.3412 


3.105 


8.17 


0.673 


2.287 


10.82 


3.94 


6.89 


0.967 


Note. - 


- The values of Pc(6) 


are in 


units of lO^g/cm 




is in units of 10^'' particles; 





7(48) and Q^^g) are in units of lO^^g • cm^; i?e(8) is the equatorial radius in 10^ cm; AMc^, 
Wr(49), and VFg(49) are in units of 10^^ erg. 
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Table 2. Isolated Magnetic White Dwarfs 



IVToTYlP 


r {pc) 


B (MC\ 

±J I iVJ-VJI I 


-'rot 


1\± 1 1\±(J\ 


PG 1031+234 


142^ 


5001^ 


3.4 hr 


0.93^ 


EUVE J0317-855 


35^^ 


450 


725 s 


1.35 


PG 1015+015 


66^ 


90 


99 min 


1.15^= 


Feige 7 


49a 


35 


2.2 hr 


0.6 


G99-47 


8^= 


25 


1 hr? 


0.7r 


KPD 0253+5052 


81^ 


17 


3.79 hr 




PG 1312+098 




10 


5.43 hr 




G217-037 


11" 


<0.2 


2-20? hr'i 


0.89 



^ Heyl (2000). 

^5 ranges from 500 to 1000. 

Liebert, Bergeron, & Holberg (2003). 
^ Downes (1986). 

^-Prot = 2 hr was used for the calculation. 



Table 3. Strain Amplitudes and Fluxes for Isolated White Dwarfs with /3 = 10 ^. 



Name /i+ Ft (Gyr) 77 



PG 1031+234 


6.0 


X 


10- 


-29 


6.1 


X 


10- 


-17 


11 


1.0 


X 


10- 


-2 


EUVE J0317-855 


1.0 


X 


10- 


-27 


6.7 


X 


10- 


-15 


1.7 


4.0 


X 


10- 


-3 


PG 1015+015 


9.3 


X 


10- 


-30 


1.1 


X 


10- 


-18 


571 


7.1 


X 


10- 


-4 


Feige 7 


1.6 


X 


10" 


-28 


4.9 


X 


10- 


-17 


125 


5.2 


X 


10- 


-4 


G99-47 


3.5 


X 


10" 


-27 


5.9 


X 


10- 


-16 


50.6 


3.7 


X 


10- 


-4 


KPD 0253+5052 


2.9 


X 


10- 


-30 


4.6 


X 


10- 


-20 


11,845 


3.5 


X 


10- 


-4 


PG 1312+098 


1.5 


X 


10- 


-30 


3.8 


X 


10- 


-21 


70,266 


2.0 


X 


10- 


-4 


G217-037 


9.0 


X 


10- 


-31 


8.2 


X 


10- 


-23 


2.4 X 10^ 


4.1 


X 


10- 


-6 



